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Driven A;-mers: Correlations in space and time 
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Steady-state properties of hard objects with exclusion interaction and a driven motion along a 
one-dimensional periodic lattice are investigated. The process is a generalization of the asymmetric 
simple exclusion process (ASEP) to particles of length k, and is called the fc-ASEP. Here, we ana- 
lyze both static and dynamic properties of the fc-ASEP. Density correlations are found to display 
interesting features, such as pronounced oscillations in both space and time, as a consequence of the 
extended length of the particles. At long times, the density autocorrelation decays exponentially in 
time, except at a special fc-dependent density when it decays as a power law. In the limit of large 
fc at a finite density of occupied sites, the appropriately scaled system reduces to a nonequilibrium 
generalization of the Tonks gas describing the motion of hard rods along a continuous line. This 
allows us to obtain in a simple way the known two-particle distribution for the Tonks gas. For large 
but finite fc, we also obtain the leading-order correction to the Tonks result. 

PACS numbers: 05.70.Ln, 05.60.Cd, 87.10.Hk 



I. INTRODUCTION 

The asymmetric simple exclusion process (ASEP), a 
paradigmatic model of nonequilibrium statistical me- 
chanics, involves hard-core particles undergoing biased 
diffusion on a lattice in the presence of an external drive 
[ll-Q. The generalization of the ASEP to an exclusion 
process of hard-core extended objects (fc-mers, each of 
which occupies k consecutive sites) is referred to as the 
fc-ASEP. It was first introduced to model protein synthe- 
sis inside living cells During the synthesis, ribo- 
somes move from codon to codon along messenger RNA, 
read off genetic information, and generate the protein 
stepwise. Modelling the codons by lattice sites and the 
ribosomes by fc-mers, we recover the fc-ASEP. The spatial 
extent of the fc-mers takes care of the blocking of several 
codons by a single ribosome; steric hindrance, which pre- 
vents overlap of ribosomes, is modelled by the exclusion 
constraint. 

Earlier studies of the fc-ASEP in one dimension in- 
volved analyzing the steady-state density profile in an 
open system [3] , the time-dependent conditional prob- 
abilities of finding the fc-mers on specific sites at a given 
time Q, the dynamical exponent Q, the phase diagram 
of the system with open boundaries [10l - [l3| . the hydro- 
dynamic limit governing the evolution of the density [l^ , 
and the effects of defect locations on the lattice on steady- 
state properties [H, [l^ . Some aspects of the fc = 2 case 
of the fc-ASEP were studied earlier in the context of a 
model of driven, reconstituting dimers [l7| . 

Here, we are concerned with the fc-ASEP on a one- 
dimensional (ID) periodic lattice. At long times, the pro- 
cess settles into a nonequilibrium steady state in which 
all configurations with a given number of fc-mers have 



equal weights In this work, our focus is on correla- 
tion functions, both static and dynamic. 

We compute static correlations in two different ways: 
(i) by counting the number of relevant configurations, 
and (ii) by mapping the fc-ASEP to a zero-range process 
(ZRP) lisll and then by employing a matrix product for- 
malism [l9[. Dynamic correlations in the fc-ASEP are 
derived by mapping the fc-ASEP to an equivalent ASEP 
with a smaller number of sites and by using the known 
dynamic properties of the latter [13 . We show that den- 
sity correlations exhibit pronounced oscillations in both 
space and time as a consequence of the extended length 
of the fc-mers. 

One may also consider the fc-ASEP in the continuum 
limit, i.e., in the joint limit of large fc and vanishing lattice 
spacing, 5 — >■ 0, while keeping the product a ~ kS and 
the density of occupied sites fixed and finite. Such a limit 
was considered previously to obtain the hydrodynamic 
behavior of the continuum system (20| . In this limit, the 
model describes hard rods of finite length a, undergoing 
biased diffusion along a continuous line in the presence 
of an external drive. In the case of unbiased motion, this 
continuum model was studied earlier by Tonks plj . This 
so-called Tonks gas has an equilibrium steady state in 
which quantities of physical interest, e.g., the equation 
of state and the two-particle distribution function, have 
been worked out exactly [111, [22| • The continuum limit of 
the fc-ASEP is a nonequilibrium generalization, and may 
be called the driven Tonks gas. 

The known two-particle distribution function of the 
Tonks gas is recovered straightforwardly by taking the 
continuum limit of the fc-ASEP, on noting that the 
steady-state measure of configurations is the same, 
whether or not the system is driven. We also obtain 
the leading-order correction to the Tonks result when fc 



is large but finite. The Tonks result, after including the 
leading-order correction in 1/fc, turns out to be a good 
approximation to the fc-ASEP equal-time spatial corre- 
lation even for not too large values of k (for example, 
fc = 13 when the density of occupied sites is 0.75). 

The paper is organized as follows. In Sec. HIl we define 
the fc-ASEP and discuss its steady-state measure. In the 
following section, we derive closed-form expressions for 
the steady-state equal-time engine-engine and density- 
density correlations. The former function describes corre- 
lation between the right end (the "engine" ) of one fc-mer 
with that of another at the same time instant. Consid- 
ering the continuum limit of the fc-ASEP, we derive the 
known two-particle distribution function for the Tonks 
gas. We also derive the leading finite- fc correction to the 
Tonks result. In Sec. IIVI we address the steady-state 
dynamics by computing the fc-mcr current and kinematic 
wave velocity associated with transport of density fluc- 
tuations. We then discuss a mapping of the fc-ASEP 
to an equivalent ASEP, and the behavior of the fc-ASEP 
temporal density-density correlation, whose scaling prop- 
erties are derived by utilizing the mapping. In Appendix 
1X1 we discuss a different method to obtain static correla- 
tions in the fc-ASEP through a mapping to an equivalent 
ZRP. 



II. THE fc-ASEP 

A. Definition 

We consider a number, N, of fc-mers that are subject to 
hard-core exclusion and are distributed on a ID periodic 
lattice of L sites. The sites are labeled by the index 
i = 1,2,...,L. Each fc-mer occupies k consecutive lattice 
sites. The fc-mers are hard objects that cannot break into 
smaller fragments. The density of occupied sites in the 
system is given by 

_Nk 

P=—- (1) 

We specify the location of a fc-mer on the lattice by the 
site index of its rightmost end, and we call this end the 
"engine" of the fc-mer. We denote the occupation of the 
ith site by rii = 1 or 0, according to whether the site is 
occupied or vacant, respectively. A fc-mer is then repre- 
sented by a string of k consecutive I's and a configuration 
of the model by an L-bit binary string composed of O's 
and I's. 

The system of fc-mers evolves according to a stochas- 
tic Markovian dynamics: in a small time dt, a fc-mer 
advances forward (respectively, backward) by one lattice 
site with probability pdt (respectively, qdt), provided the 
site is unoccupied. The dynamics conserves the total 
number of fc-mers in the system. The elementary dy- 
namical moves may be represented as 
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pdt 

(111. ..1)0 ^ 0(111. ..1), 

qdt 
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FIG. 1: (Color online) The fc-ASEP on a ring, showing the 
allowed and disallowed dynamical moves of trimers (fc = 3). 
Here, fc-mers are represented as fc connected circles. The 
rightmost end of a fc-mer (the "engine") is indicated by a 
filled red circle. 



where the fc-mcr has been represented by a string of fc 
consecutive I's enclosed by brackets. Figure [1] illustrates 
the allowed and disallowed moves of trimers (fc = 3). 

For p ^ q, the fc-mers move preferentially in one di- 
rection along the lattice, and the system at long times 
reaches a noncquilibrium steady state with a steady cur- 
rent of fc-mers. 

When p = q, the fc-mers diffuse symmetrically to the 
left and to the right. In this case, the model is a general- 
ization of the symmetric simple exclusion process (SEP) 
of hard-core particles to fc-mers, and may be referred to 
as the fc-SEP. At long times, the fc-SEP settles into an 
equilibrium steady state. 



B. Relation to a model of diffusing, reconstituting 
dimers 



Earlier studies of a variant of the fc-ASEP dealt with 
dimers (fc = 2) that, in general, do not retain their identi- 
ties and are allowed to reconstitute [l3, [111 . Specifically, 
on a lattice of L sites, single particles or monomers (de- 
noted by 1) and paired particles or dimers (denoted by 
11) are distributed with at most one particle per site. 
The dynamics in a small time dt involves a dimer mov- 
ing by one lattice site, either forward with probability pdt 
or backward with probability qdt, without violating the 
hard-core constraint on site occupancies. The pairing of 
the dimers is impermanent, thereby allowing for recon- 
stitution. For example, in the sequence of transitions, 
11010 OHIO 01011, the middle particle is paired 
with the particle to the left in the first transition and 
with the particle to the right in the second transition. 

Both in the symmetric {p — q) [23| and in the asym- 
metric {p 7^ q) [Ol case, the phase space of the sys- 
tem breaks up into an infinite number of dynamically 
disjoint sectors. A non-local construct, called the ir- 
reducible string (IS), uniquely labels the different sec- 
tors. The IS for a given configuration is constructed from 
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the corresponding L-bit binary string by deleting recur- 
sively any pair of adjacent I's until no further deletion 
is possible. The model exhibits dynamical diversity with 
quantities like the density autocorrelation showing strong 
sector-dependent behaviors that range from power laws 
to stretched exponentials. 

The problem of hard non-reconstituting dimers, i.e. 
the case fc = 2 of the fc-mer system considered in this 
work, corresponds to a particular sector of the reconsti- 
tuting dimer problem, namely, the one with the null IS 
0000 . . . 000. It is easy to check that this is the sector in 
which dimers do not reconstitute. 



site is occupied by an engine while the following site is 
vacant; this is given by the ratio L-(fe+i)/^- 
thermodynamic limit, i.e. in the limit N — > oo,L — > oo, 
while keeping k and p fixed and finite, this joint proba- 
bility equals p{l — p)/[p + k{l — p)]. 

We note that for the /c-SEP, the condition of detailed 
balance implies that all configurations with a given num- 
ber of /c-mers have equal weights in its equilibrium steady 
state. It thus follows that the exclusion process with k- 
mcrs has the same steady state, irrespective of whether 
the fc-mers have a biased or an unbiased motion. 



C. The steady state 

In the steady state of the fc-ASEP, every microscopic 
configuration C with a given number of fc-mers is equally 
likely, and hence, occurs with probability l/il, where is 
the total number of configurations [ll| . In such a state, 
for every transition away from C to another configuration 
C", it is simple to construct a distinct and unique config- 
uration C" that evolves to C at the same rate, thereby 
ensuring stationarity. This argument holds for periodic 
boundary conditions and fails in the case of open bound- 
aries where the exact steady state is hard to obtain and 
is as yet unknown. 

Now, n is determined by counting the different possible 
ways of distributing a number, N, of fc-mers over a lattice 
of L sites with periodic boundaries. First consider free 
boundaries, in which case we have an open chain of L 
sites on which the fc-mers are placed. Then, the number 
of possible ways of distributing the fc-mers is 



n 



N.L 



L- Nk + N 
N 



(3) 



For a periodic lattice, the number of ways is determined 
by first considering all configurations in which an arbi- 
trary but fixed site remains occupied by an arbitrarily 
chosen but fixed engine. The number of such configura- 
tions is simply However, the engine could be 
chosen to be one of the N identical engines available and 
may be placed over any of the available L sites. Thus, 
the total number of distinct configurations is [ll| 



III. STEADY STATE STATICS 

A. Engine-engine correlation 

Let Ei denote the engine occupation variable for site 
i, taking values I or according to whether the site is 
occupied by an engine or not, respectively. The (un- 
subtracted) equal-time engine-engine correlation func- 
tion Eu {r) is defined as 



£k{r) = {E,E, 



i+r I 



(5) 



where the angular brackets denote averaging with respect 
to the steady state. Note that £k{r) is identically zero 
for r < k. 

Now, ffc(r) for r > k has non-zero contributions from 
all configurations in which there are two engines at the 
i-th and {i + r)-th sites, with the gap between the two 
engines containing r — k sites that are occupied by any 
number, m, of fc-mers between and the maximum num- 
ber that may be placed in the gap, while the left over 
N — (m -|- 2) number of fc-mcrs is distributed over the 
remaining L — (r — k + 2k) sites. The maximum number 
of fc-mers that may be placed over the gap of r — fc sites is 
given by [(r — fc)/fcj , where [x\ denotes the floor function 
that gives the largest integer not greater than x. Noting 
that in the steady state of the fc-ASEP, all configurations 
have equal weights, we get 



L(r-fe)/feJ 



Qlrcc Qt] 



(6) 



-n 



N 
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N-l,L-k 
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Nk + N 
N -1 



(4) 



The steady-state probability that a randomly chosen 
site is occupied by an engine is obtained by considering 
all possible ways of distributing a number, iV — 1, of fc- 
mcrs over a lattice oi L — k sites with free boundaries. 
Since all configurations of the fc-ASEP are equally likely 
in the steady state, the desired probability is given by the 
ratio ^^Ipli L-k/^ ~ P/^- Similarly, the probability that 
a randomly chosen site is vacant in the steady state may 
be shown to be equal to I — p. Next, we consider the joint 
probability that in the steady state, a randomly chosen 



In the thermodynamic limit, one can show by using 
Eq. (jSj that, for arbitrary integers mi and TO2, 



n 



free 



Ofree " ^0 U - Po) ) [' ) 



where 



Po 



k{l-p) + p 



(8) 



Later, in Sec. IIVI by using a mapping of the fc-ASEP to 
the usual ASEP, we show that the quantity po is in fact 
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FIG. 2: (Color online) The fc-ASEP density-density correla- 
tion Cfe(r) and engine-engine correlation Ekij'), evaluated nu- 
merically by using Eqs. (|12[l and ((9| for fc = 64 and p = 0.75. 
The inset shows that Vk{r) and - pf j p^\Ek{r -|- fc - 1) 

are equal, in agreement with Eq. (|17p. 



the particle density in the latter. Using Eq. (O, we find 
that in the thermodynamic limit, Eq ([6]) reduces to 



£M = 



L(r-fc)/fcJ 



pI 



E 



l + po(fc-l) 
— fc — km 4 



— k—kn 



> fc. 



(9) 



For fc = 2, the sum in Eq. ([9]) can be evaluated exactly 



to obtain £k{r) 



Po 



[1 - (-Po)''-'] [i3- For fc > 3, 



the function £k{r) may be numerically evaluated by us- 
ing Eq. dni). Figure [5] shows the result for fc = 64. We see 
that the engine-engine correlation exhibits damped oscil- 
lations in space, a hallmark of systems with hard-core 
interactions between the constituents 2J], where purely 
entropic considerations apply. For example, the fact that 
£k{r) has its first minimum at r = 2fc — 1 and rises again 
at r = 2k may be understood to be due to the possibil- 
ity that when r = 2k, there may be configurations with 
an additional engine between the two engines that are 
occupying sites i and i + 2k. 



B. Density-density correlation 

The (unsubtracted) equal-time density-density corre- 
lation function Ck{r) is defined as 



Cfc(r) = {niHi+r)- 



(10) 



Here, rij is the occupation variable for site i, taking values 
1 or 0, according to whether the site is occupied or is 
vacant, respectively. Evidently, Ui can be expressed in 



terms of the engine occupation variable Ei as 

fe-i 

rii = ^ Ei+rn- (11) 
?n— 



Using the above equation, Ck{r) may be expressed in 
terms of the equal-time engine-engine correlation as 



fc-i 

Cfc(r) = k£k{r)+Y^ m\£k{r + k-m)+£k{r-k + m) 

m— 1 

(12) 

where it is understood that £k{r) is zero for r < k. 

Alternatively, Ck {r) can be computed in a straightfor- 
ward way from the equal-time vacancy-vacancy correla- 
tion, defined as 



Vkir) = (riin^+r), 
where rii ~ I — ni , so that 

Ck{r) = Vk{r) + {2p-\). 



(13) 



(14) 



Now, Ffc(r) for r > 1 has non-zero contributions from 
all configurations in which both the i-th and {i -\- r)-th 
sites are vacant, with the gap between the two containing 
r — 1 sites that are occupied by any number, m, of fc-mers 
between and the maximum number that may be placed 
in the gap, while the left over N number of fc-mers is 
distributed over the remaining L — (r — 1 + 2) sites. We 
get 



[(r l)/k\ Qfrcc Qfree 

Vkir) = 2^ 



m=0 



(15) 



In the thermodynamic limit, we find that 
(1-po)^ 



Vkir) = 

L(r-l)/fcJ 

E 

m=0 



1 + Poik-l) 
r — 1 — km + m 



PS'i^-Po)'^ 



r>l. 



(16) 



It follows from the definition that 14(0) = 1 — p. On 
comparing Eq. p6)) with Eq. ([9]), and noting that (1 — 
Po)Vpo = ^'(1-p)Vp', we get 

Vkir)= ^ , ^' £kir + k~l); r > 1. (17) 
P 

Using Eq. (fT4| . we get 

Ckir) = ^'^^r^^' gfc(r + fc-l) + (2p-l); r > 1. (18) 
P^ 

For r = 0, we have C(0) = p. 

Figure [5] shows Ckir), computed using Eqs. (jl2p and 
dni), for fc = 64 and p = 0.75. The oscillations in Cfe(r) 
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FIG. 3: (Color online) Scaling approach of the engine-engine 
correlation £k(r) to the Tonks limit; k'^£k{r) vs. r/k at fixed 
p = 0.75 and a — 400 shows data collapse for large k according 
to Eq. p3ll . 




FIG. 4: (Color online) Comparison of k'^£k{r) with the Tonks 
result, a^g{r/k), for k = 13, a = 400, and p = 0.75, show- 
ing the discrepancy between the two for large but finite k. 
The discrepancy is resolved on including the correction to the 
Tonks result, to leading order in 1/fc, as is shown by a com- 
parison of k^£k{r) with the function a^g(r/k) + {p/k)h(r /k) 
(Eq. (123)). 



may be related to those in £k{r) by using Eq. p^ . For 
instance, since £fe(r) has its first minimum at r = 2fc — 1, 
it follows that the first minimum of Cfc(r) occurs at r = 
k. The inset of Fig. H shows that Vk{r) and [^^(1 - 
p)'^ / p'^]£k{r + fc — l)forr>l are equal in accordance 
with Eq. (fTT)) . 

Later, in Appendix [Al we discuss an alternative 
method to obtain static correlations in the fc-ASEP 
through a mapping to a zero-range process and employ- 
ing a matrix product formalism [l9[. 



Then, if a is the rod length, it is known that [12] 

P{R) = g{x), 
q{x) = — 7 > Aix — mj- -— — 

xexp(--^), (19) 



C. Continuum limit: Driven Tonks gas 



In a suitable continuum limit of the fc-ASEP, discussed 
below, the model reduces to one of hard rods, which have 
exclusion interaction, and which are undergoing driven, 
diffusive motion along a continuous line. In the absence 
of drive, this continuum model was studied by Tonks as 
a ID interacting system with an equilibrium steady state 
in which thermodynamic properties like the equation of 
state can be worked out exactly [2l|. When the motion 
is driven, the continuum limit of the fc-ASEP becomes 
the driven Tonks gas. 



A quantity of physical interest for the equilibrium 
Tonks gas is the two-particle distribution P(i?i,i?2), de- 
fined such that P{Ri, R2)dRidR2 is the joint probabil- 
ity of finding the rightmost end of one hard rod be- 
tween i?i and Ri -|~ dRi and that of another between 
R2 and i?2 + dR2 . For a translationally invariant system, 
P(i?i, i?2) is a function of the separation R= \R2 — Ri\. 



where x = R/a is a reduced distance. Here, I = l/apx, 
where px is the density of rods, and A{x) is the unit step 
function: 



Aix) 



for a; < 0, 

1 for a; > 0. 



(20) 



We now show that the continuum limit of the fc-ASEP 
engine-engine correlation easily yields Eq. (jl9p . Such 
a derivation is justified by the fact that, as discussed in 
Sec. Ill CI the fc-ASEP has the same steady-state measure 
of configurations for both unbiased and biased motion of 
the k-mers. This fact further implies that Eq. ([T^ also 
holds for the driven Tonks gas. 

The continuum limit of the fc-ASEP is obtained by 
considering the joint limit fc — > 00, r — > 00, and the 
lattice spacing (5 — > 0, while keeping R = rS, a = k6, and 
p fixed and finite [l^l • The fc-ASEP then describes biased 
motion of hard rods of length a along a continuous line. 
The density of hard rods is pT = p/a. In this limit, when 

Po = pS/[a{l - p)] and (1 - po) = exp - pS/[a(l - p)] 
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Eq. ^ reduces to 

''-^ ' la^ ^ m!(Z- l)™+i 

m=0 ^ ' 

( i?/a-(m+l)\ 

xcxp(^ r~i — /■ ^ ^ 

Comparing the right-hand side of the last equation with 
Eq. ([T^ . and noting that Rja — r jk, we find that 

'52.9(9 • (22) 

Now, since (5 = a/fc, we find that in the continuum hmit, 
i.e., in the hmit fc— >-oo,r— >oo,(5— >0, while keeping 
a, i?, and p fixed and finite, Ekif) for different fc has the 
scaling form 

^.(0 = ^,9©. (23) 
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FIG. 5: (Color online) Scaling approach of the density-density 
correlation Cfc(r) to the Tonks limit: Ck{r) vs. r/k -I- 1 at 
fixed p = 0.75 and a = 400 exhibits data collapse for large 
k in accordance with Eq. p4|l . The data are obtained by 
numerically evaluating Eqs. ([5} and (|12|l . 



Moreover, in the continuum limit, defining P{R) = 
£k(r)/5'^, we find from Eq. ((22)) that is precisely 

in the form of Eq. ([T^ . We have thus derived the two- 
particle distribution, valid for both the equilibrium and 
the driven Tonks gas, by considering the continuum limit 
of the fc-ASEP. 

Figure [3] shows plots of k'^£k{r) for different k at fixed 
p and a, evaluated using Eq. ([5]). As fc increases, the 
curves show a good scaling collapse, in accordance with 
Eq. ([23)) . For large but finite fc, the right-hand side 
of Eq. ([231) has finite-fc corrections, as is suggested by 
the discrepancy between k^£k{,f) and a^g{r/k), shown in 
Fig. 131 The leading-order correction to Eq. (P3| will be 
discussed in the following subsection. 

The behavior of the density correlation Cfe(r) in the 
continuum limit may be easily obtained by using Eqs. 
(IT51) and (1121). We find that in this hmit, for different k 
at fixed p and a, we have 



-9 



(^ + l)+(2p-l). 



(24) 



Figure |S] shows plots of Ck{r) for different fc at fixed p 
and a, evaluated using Eqs. ^ and (|T2l) . As k increases, 
the curves show a good scaling collapse, in accordance 
with Eq. 



Finite-fc corrections to Tonks two-particle 
distribution 



D. 



In order to compute finite-fc corrections to Eq. 
we evaluate the engine-engine correlation for finite r, 
fc, and (5, with r ::g> 1, fc ^ 1, and (5^1, keep- 
ing R = rS,a = kS, and p fixed and finite. Then, 



on substituting po 



and (1 - po) 



exp 



ing terms to leading order in 1/fc, we get 



(^ 2k'^(i-nVi ) i'^to Elq. (0), and keep- 
order 



where 
/i(x) = 



^(a: — to) 

m— 1 

rTO(TO — 1) 



(a; - to)™+i 



exp 



(m - 1)!(/ - 1, 

TO X — TO 



/ - 1 



(26) 



L2(a;-m) {I - 1) 2{l - I f 



Here, A(x) is the unit step function defined in Eq. (|20p . 
Figure m shows that inclusion of the leading-order correc- 
tion to the Tonks result, as in Eq. (pS)) . indeed resolves 
the discrepancy between k'^Ekir) and o?g{r /k). 



IV. STEADY STATE DYNAMICS 

A. Current and kinematic wave velocity 

In discussing the current in the system, we need to 
distinguish between that associated with the motion of 
engines, and that with the fc-mers. Contributions to the 
engine current across a bond (i, arise when either (i) 
the i-th site is occupied by an engine, while the {i + l)-th 
site is vacant, or, (ii) the (i -I- l)-th site is occupied by an 
engine, while the (i — fc+l)-th site is vacant. On using the 
results of Sec. Ill CI we find that in the thermodynamic 
limit, the average engine current in the steady state is 
given by 



Je 



{p-q)p{l - p) 
p + k{l-p) 



(27) 
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To compute the fc-mer current. J, note that associated 
with the motion of the engine to an adjacent site is the 
shding of the corresponding fc-mer across (fc— 1) bonds, so 
that J = fcjg. It can be checked that J has a maximum 
at the density = Vk/^^/k + 1). 

The kinematic wave velocity vk = dJ / dp accounts for 
the transport of density fluctuations through the system 
in the steady state [25] . Wc find 



VK = k(j) - q) 



(fc-l)p2 + fc(l_2p) 



[p + k[l-p)f 



(28) 



Evidently, vk vanishes if the density is pc- As wc dis- 
cuss below, Vk plays an important role in determining 
the form of the temporal decay of the density autocorre- 
lation. 



B. Mapping to the ASEP: Wheeling velocity 



fc-ASEP configuration ASEP configuration 
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ASEP site that moves around the ASEP ring with a finite 
mean velocity. This phenomenon is known as wheeling, 
and the mean velocity is called the wheeling velocity W 
[l7| . This velocity plays an important role in the scaling 
properties of the temporal density-density correlation of 
the fc-ASEP, as we discuss in the next subsection. 

The mapping involves representing each fc-mer by a 
hard-core particle that corresponds to the engine of the 
fc-mer, as illustrated in Fig. [Sj In this way, every fc-ASEP 
configuration is mapped to a unique configuration in the 
ASEP. Associated with the motion of a fc-mer is that of 
the corresponding ASEP particle according to the ASEP 
dynamics. We now see that the quantity po in Eq. (jH)) is 
the particle density in the equivalent ASEP. 

It can be seen from Fig. [6] that in the process of map- 
ping, the ASEP image of a fixed fc-mer site moves around 
the ASEP lattice as a result of the fc-mer motion. For ex- 
ample, consider the transition (111 . . . 1)0 — !• 0(111 ... 1). 
It is easy to sec that in this transition, the ASEP images 
of fc-ASEP sites containing the and the leftmost 1 in 
the string (111. ..1)0 do not change, while the images 
of those containing the remaining I's in the string in- 
crease by one unit. On the other hand, in the transition 
0(111 . . . 1) ^ (111 . . . 1)0, the ASEP images of fc-ASEP 
sites containing the and the rightmost 1 in the string 
0(111 ... 1) do not change, while the images of those con- 
taining the remaining I's in the string decrease by one 
unit. This motion of an ASEP site corresponding to a 
fixed site in the fc-ASEP is the phenomenon of wheeling. 
As a result, the displacement of the ASEP image of a 
fixed fc-ASEP site in time t is given by 



Ar(t) = Wt + (j){t), 



(29) 



where 0(t) is a random variable with zero mean, arising 
from the stochasticity in the dynamics. Referring to the 
results on joint occupation probabilities in Sec. Ill C[ we 
find that the wheeling velocity is given by 



W 



(p-g)(fc~l)p(l-p) 
p + fc(l-p) 



(30) 



C. The temporal density-density correlation 

The temporal density-density correlation function in 
the steady state of the fc-ASEP is defined as 



FIG. 6: (Color online) Mapping of a typical fc-mer configura- 
tion to an ASEP configuration and their subsequent evolution 
in time. The first site of the ASEP lattice may be defined in 
more than one way; the figure illustrates one possibility. 

We now discuss a mapping of the fc-ASEP with a num- 
ber, N , of fc-mers on a ID periodic lattice of L sites to an 
ASEP of N hard-core particles on a ID periodic lattice of 
L' = L — N{k — 1) sites. We show that as a result of the 
mapping, a fixed site in the fc-ASEP corresponds to an 



Cfc(r= |z-j|,t) = (n,(0)n,(t))-p2 



(31) 



where ni(t) denotes the occupation index of site i at time 
t. In particular, the density autocorrelation is given by 
Ckit)=Ck{0,t). 

We study Ck (t) by performing Monte Carlo simulations 
of the fc-ASEP in the steady state. Figure [7] shows that 
Ck{t) oscillates in time. The inset shows that for large fc, 
the autocorrelation for different fc at a fixed density p is 
initially a function oit/k. This behavior holds up to a fc- 
dependent time. To understand the dependence of Ck{t) 
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0.16 




Time t 



FIG. 7: (Color online) Short-time behavior of the fc-ASEP 
density autocorrelation Ck{t) for fc = 60 and p = 0.8. The 
data are obtained from Monte Carlo simulations. The inset 
shows Ck{t) vs. t/k at fixed p = 0.8, illustrating scaling for 
large k. 



1 




t 

FIG. 8: (Color online) Long-time decay of the fc-ASEP density 
autocorrelation Ckif), as a power law in time at the compen- 
sating density poc, and as an exponential at other densities. 
The data are obtained from Monte Carlo simulations. Here 
fc = 4, so that poc = 1/3. The black line has the slope of 
-2/3. 



on the ratio t/fc, wc note that at short times, the relevant 
time scale is set by the time T(fc) that an occupied site 
takes to fall vacant. This is consistent with Ck{t) being 
a function of t/rik). An upper bound on T(fc) may be 
estimated as the time Te(fc) « fc/we, the time that an 
engine takes to move by its own length. Here Ve is the 
velocity of an engine, which may be obtained from Eq. 
(f27| as Ve = ''p+fc(il~f ^ ■ In the limit of large fc, one 
finds that re(fc) sa fc, so that Ck{t) is a function of t/k, as 
observed. 

We now discuss the behavior of Ck{t) at long times. 
To proceed, we examine the function Ckir^t) for which 
an earlier study in the case k = 2 has illustrated that 
in the limit of long times and large distances, it assumes 
a particular scaling form To obtain the scaling for 
general k, we utilize the mapping to the ASEP discussed 
in Sec. IIVBI and invoke known scaling properties of the 
density correlation in the latter. 

In the ASEP (fc ~ 1), the temporal density-density 
correlation function, Ci(r = |i — = {ni{0)nj{t)) — /Oq, 
in the scaling limit follows the form [2^ 



Ciir,t) (X t-^/^F{u); 



(32) 

Here, Jq and vqk are, respectively, the steady-state cur- 
rent and the kinematic wave velocity in the ASEP, given 

by Jo = (p-g)Po(l-po), andi;oif = (p-g)(l-2po) H- In 
the limit of large i t, it is known that F{u) ~ exp(— ^|up), 
with -0.295 [2^. 

It is evident from the fc-ASEP to ASEP mapping dis- 
cussed above that at long times, neglecting the stochas- 
tic part (j){t) in the displacement of a mapped site in the 
ASEP, the correlation Ck{r,t) has a behavior similar to 
Ci (r + Wt, t) [13] ■ Thus, for a fixed k, in the limit of long 



times and large distances, Cfc(r, t) follows the scaling form 

Ck{r,t) oc t-^/^F{uy, u' = l(J„t^yy^(^r+iW-voK)t 

(33) 

The autocorrelation Ck{t) behaves asymptotically as 



(34) 



where k is a constant determined by the difference 
{W — Vqk)- Thus, at long times, Ck{t) decays as an ex- 
ponential in time, unless the density p is such that the 
difference vanishes. In this case, the autocorrelation at 
late times decays in time as a power law: Ck(t) ~ t~'^^^. 
The corresponding ASEP density is called the compen- 
sating density poc [l7j . and satisfies 



pl,{k - 1) + 2poc -1 = 0. 
Solving for the positive root, we get 

1 



PQc 



Vk + l' 



(35) 



(36) 
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which matches with the result for k — 2 derived in 
Note that corresponding to poc is the fc-ASEP density of 
occupied sites Pc, mentioned in Sec. IIV Al at which the 
fc-mer current J is maximized and the kinematic wave 
velocity vk is zero. 

Figure H] shows Ck {t) as a function of time for three 
values of the ASEP density, namely, the compensating 
density poc and two other values on either side. We see 
an asymptotic t~^^^ decay of the autocorrelation at the 
compensating density and an exponential decay at other 
densities, in accordance with our analysis above. 
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FIG. 9: (Color online) Mapping of the fc-ASEP to the ZRP. 
The fc-ASEP vacancies are considered as sites in the ZRP, 
while an uninterrupted sequence of fc-mers in front of a va- 
cancy is regarded as a set of particles occupying the corre- 
sponding ZRP site, such that the number of particles equals 
the number of fc-mers in the sequence. 



The delta function stands for overall particle conserva- 
tion. 

The steady state weight of any fc-ASEP configuration is 
obtained by mapping it to a configuration in the equiva- 
lent ZRP and computing its weight by utilizing Eq. (jAip . 
In our case, since /(m) = 1 for all m, all /c-ASEP config- 
urations with a given number, N, of fc-mers are equally 
likely. Now, using the formalism discussed in [l9| . one 
may rewrite the steady-state weight of any fc-ASEP con- 
figuration (ni,n2, . . . ,nL) in a matrix product form by 
replacing each occupation number rii by either a matrix 
_D or a matrix E depending on whether rii is 1 or 0, 
respectively. From the correspondence between the fc- 
ASEP and the ZRP, we have 
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Appendix A: Mapping of the fc-ASEP to the ZRP 

In this appendix, we discuss a mapping of the fc-ASEP 
to a zero-range process (ZRP) and show how static corre- 
lations in the former are obtained by using the mapping. 
This method of obtaining the correlations is an alterna- 
tive to the one discussed in Sec. IIIII 

In the ZRP, an unrestricted number of particles resides 
on lattice sites and hops between sites with a rate that 
depends only on the number of particles on the departure 
site Each fc--ASEP configuration can be mapped to 
a unique ZRP configuration in the following way: one 
considers vacancies (O's) in the fc-ASEP as sites in the 
ZRP and an uninterrupted sequence of fc-mers following 
a vacancy as particles residing on the corresponding ZRP 
site, with the number of particles equal to the number of 
fc-mers in the sequence. The ZRP has N particles and 
M = L — Nk sites labeled by the index i. A generic ZRP 
configuration is the set {rui} = {mi, 1122 . . .uim)-, where 
rrii is the number of particles on the i-th site (see Fig. 
[9]). The motion of a fc-mer in the fc-ASEP translates to 
hopping of a particle from a ZRP site to its right or left 
neighbor with rates p or respectively. 

For an arbitrary hop rate, the ZRP has a product mea- 
sure stationary state In our case, where the rates do 
not depend on the number of particles on the departure 
site, the steady-state weight of any configuration {rrii} is 

M 



P({mO) = Tv[ED 



km I 



.ED 



km,} 



']s{ 



N 



M 



,(A2) 



where Tr denotes the usual matrix trace operation. 

Without loss of generality, one may take E = \a){/3\, 
where the vectors \a) and (/3| are to be determined. This 
choice of E together with Eqs. (jA2p and (jAip demand 
that for any positive integer m, the matrix D satisfies 



\D 



mk 



a 



f{m) 



1. 



(A3) 



Also, {P\D^\a) = for positive integers j which are not 
multiples of fc. A simple fc-dimensional representation of 
matrices E and D is when they have non-zero elements 
Ell = 1 and D^i = 1 = Di i^i, that is. 



E = \a){/3\, with \a) ^ 



(II 



fe-i 



(A4) 



Here, the set represents the standard basis vectors in 
fc-dimensions. This choice ensures that the weight of any 
configuration with one or more blocks of / particles is zero 
if / is not an integral multiple of fc; all other configurations 
are equally probable. 

Let us mention that the matrix formulation discussed 
above is different from the Matrix Product Ansatz 
(MPA) of Derrida et al. in which matrices satisfy 

specific algebraic relations dictated by the system dy- 
namics. For models with ZRP correspondence, the ma- 
trices generically satisfy Eq. (|A3|) . and therefore, depend 
only on the weights f{m). It is always possible to get one 
representation of these matrices, whereas finding explicit 
representation of the MPA matrices is non-trivial. 

Partition function Zl(z): The first task in comput- 
ing fc-ASEP static correlations is to find the partition 
function of the system, which is conveniently done in 
the grand canonical ensemble by associating the fugac- 
ity z with any occurrence of the matrix D. This gives 
Zl{z) = Tr[C^], where C = zD + E. The configuration 
with no vacant site is not dynamically accessible. Thus, 
Zl{z) is given by weights of all configurations with at 
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least one vacant site: 



Zl{z) = ^Tr [(zL>)"-i£;C^-"] 

n=l 
L 

= Y^{P\C'^-{zDr''\a). (A5) 



n=l 



To proceed, we use the following generating function; 

oo 



L = l 



'I-7CI-7ZD' 



7[l + (fc-l)(7z)1 
[l-7-(7z)'^-][l-(7z)'^-] 



(A6) 



where I is the fc-dimensional identity matrix. Note that 
Z{z,^) may be interpreted as the partition function in 
the variable length ensemble. The parameters z and 7 
together determine macroscopic observables like the den- 
sity of occupied sites and the average system size. The 
density of occupied sites is p = 1 — {fii), where 



7 



1 



1 - (7z)* 



'I-7C" 



l + (fc-l)(7z)'=' 



(A7) 



The average system size is given by {L) = ^ 7^ 



z 97 • 



One may check that {L) has the radius of convergence 
. 1 



7 



(A8) 



Thus, the thermodynamic limit (L) — > 00 is achieved at 
z = z* , when every observable of the system becomes a 
function of 7 only. For example, p is obtained from Eq. 
7| as p = 7+^^(7-7) ' '^hich may be inverted to obtain 

k{l-p) 



7 



p + k{\-p) 



1 - Pe- 



rn 



Equal-time correlations Vk(r),£k{f): One may com- 
pute Vk{r) by writing it in terms of matrices as follows: 



Vk{r) = 



7 



r+l 



]C^-'\a){(3\-, 



1 



Y{l-pmC^-'\a}. (AlO) 



Now, for any integer j, one may check that 



m=0 



z™", (All) 

m / 



which results in 



L(r-l)/fcJ 
m=0 

The equal-time engine-engine correlation £k (r) may be 
similarly calculated by using the matrix formulation: 



7 



r+k 



L(r-fc)/fcJ 



L(r-fe)/feT=° 



In the thermodynamic limit, using Eq. (|A8|) . we get 

r+i L('-i)/fcJ 



7 



7 + k{l — 7) 



E(-)frcc 
")ri,r-l 



m=0 



1-7 



7 + A:(l - 7) 



m=0 



7 



Replacing 7 by 1 — po in the expressions on the right 
reduce them to those in Eqs. ([T5| and respectively. 
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